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ABSTRACT 

We demonstrate that it is crucial to account for the evolution of the starburst 
population in order to derive reliable numbers of O stars from integrated spectra for 
burst ages t > 2 — 3 Myr. In these cases the method of Vacca & Conti (1992) and Vacca 
(1994) systematically underestimates the number of O stars. Therefore the current 
WR/0 number ratios in Wolf-Rayet (WR) galaxies are overestimated. This questions 
recent claims about flat IMF slopes {a ~ 1-2) in these objects. If the evolution of the 
burst is properly treated we find that the observations are indeed compatible with a 
Salpeter IMF, in agreement with earlier studies. 

Including recent predictions from non-LTE, line blanketed model atmospheres which 
account for stellar winds, we synthesize the nebular and WR He II A 4686 emission in 
young starbursts. For metallicities 1/5 Zq < Z < Zq we predict a strong nebular B.e II 
emission due to a significant fraction of WC stars in early WR phases of the burst. 
For other metallicities broad WR emission will always dominate the He II emission. 
Our predictions of the nebular He II intensity agree well with the observations in WR 
galaxies and an important fraction of the giant H II regions where nebular He II is 
detected. We propose further observational tests of our result. 

Subject headings: galaxies: starburst — H II regions — stars: Wolf-Rayet 
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1. IMF determinations in young starbursts 

To derive the initial mass function in starburst 
galaxies is one of the fundamental goals in stud- 
ies of starforming regions. Due to their large 
luminosity, massive stars are directly visible in 
the integrated spectrum of young starbursts and 
thereby provide a unique opportunity to study 
their stellar content. In a subset of emission line 
galaxies - often referred to as Wolf-Rayet galax- 
ies (cf. Conti 1991) - broad stellar emission lines 
in the optical (most prominently He II A 4686) 
testify to an important population of Wolf-Rayet 
(WR) stars, which are the descendents of the 
most massive O stars (see e.g. Maeder & Conti 
1994). The nebular emission lines, on the other 
hand, are due to the exciting stars present in the 
starburst including the less massive OB stars. 
Observations of WR features and nebular lines 
thus contain information about stars from differ- 
ent mass ranges. WR galaxies therefore in par- 
ticular offer a unique opportunity to probe the 
upper part of the IMF in young starbursts. 

1.1. Difficulties with the non-evolving clus- 
ter approach 

Since the first discoveries of WR stars in emis- 
sion line galaxies by Allen, Wright, & Gross 
(1976) and Kunth & Sargent (1981), the recent 
work of Vacca & Conti (1992, hereafter VC92) 
provides the most detailed quantitative study 
of the massive star population in WR galax- 
ies. Their technique, developed in more detail by 
Vacca (1994, hereafter V94), allows to determine 
the relative number of WR and O stars. Com- 
parisons with stellar evolution models show that 
the large WR/0 ratios can only be explained 
by bursts of star formation occurring over short 
timescales compared to the lifetime of massive 
stars (Arnault, Kunth, & Schild 1989, VC92, 
Meynet 1995). 

At first sight the comparison of the WR/0 ra- 
tios of VC92 and Contini, Davoust, & Considere 
(1995, hereafter CDC95) with the evolutionary 



calculations of Meynet (1995), however, seems to 
show that a flat IMF slope {a ^ 1-2, compared to 
the Salpeter value of a = 2.35) is required to ex- 
plain the observations (Meynet 1995, 00095)0. 
This result is, however, in disagreement with sev- 
eral other methods, such as direct stellar counts 
in clusters (a ~ 2 - 3.1; cf. Maeder & Conti) and 
population synthesis techniques (Olofsson 1989, 
Mas-Hesse & Kunth 1991, Kriiger et al. 1992, 
Vacca et al. 1995), which show that observa- 
tions of similar (and partly even the same) ob- 
jects are compatible with a Salpeter IMF over a 
large range of metallicities. 

It can, indeed, easily be traced back why ap- 
plying the method of VC92 and V94 can lead 
to erroneously large WNL/0 ratios, which would 
imply a flat IMF. From the ratio of the luminos- 
ity in the WR feature to the H/3 luminosity, they 
derive N-sj^nl/Nqjy, the number of WNL stars 
over the number of so-called "equivalent 07V 
stars" , which are required to produce the Lyman 
continuum flux measured from the H/3 line. They 
then deduce the ratio of WNL over OV stars from 
AVnl/A^OV = ^^'wNL/A^oyv) where tjq is the 
IMF averaged ionizing Lyman continuum lumi- 
nosity of a population normalized to the output 
of one "equivalent" 07V star (see V94). 

The fundamental assumption entering the deri- 
vation of A^wnlZ-^oV) or equivalently rjo, is that 
the Lyman continuum luminosity is produced by 
an unevolved Zero-Age Main Sequence (ZAMS) 
population. Obviously this assumption is not 
compatible with the presence of WR stars — 
the descendents of O stars — which clearly indi- 
cate ages of typically at least 3 Myr. A^wnl/A'^ov 
therefore denotes the number ratio of WNL stars 
with respect to a purely theoretical ZAMS OV 

■*It should be noted that in Fig. 3 of Maeder & Conti (1994) 
the slope of the IMF has inadvertently been mislabeled, 
due to a confusion in the notation. Their values of F = — 1 , 
and —2 should correctly refer to a = 1 and 2 in the above 
notation, where the Salpeter slope is given by a = 2.35. 
See also the original paper by Meynet (1995) who uses 
a = X. 
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star population, which: 1) has no observational 
correspondence, and 2) cannot directly be com- 
pared to WR/0 number ratios derived from evo- 
lutionary models. 

In an instantaneous burst, due to the evo- 
lution of the O star population and the subse- 
quent disappearance of the most massive stars, 
the number of Lyman continuum photons, Qq, 
strongly decreases after 2-3 Myr. This is illus- 
trated in Fig. H, where we plot the time evolution 
of 

_ Cc^{M)Q,{M,t)dM 

with respect to its Z AMS value rj^ {t = 0) in an in- 
stantaneous burst. Following V94, <I> is the IMF 
with the mass limits Miow and Mup, QQ{M,t) 
is the ionizing luminosity as a function of initial 
mass and time, and logQg*^^ = 49.05. The inte- 
gral in the denominator counts the number of O 
stars defined as having T^s > 33000 K. 

The strong decrease of i]o after ~ 2 Myr im- 
plies that the required number of O stars will 
be systematically underestimated (and hence the 
WR/0 ratio overestimated) if one does not ac- 
count for the evolution of the population and uses 
the ZAMS values of a "non-evolving cluster" . For 
younger ages the approach of V94 should, how- 
ever, be valid. At the low metallicities typical 
of the objects from VC92, the decrease of r]o 
in the WR rich phases (see Fig. [l|) implies that 
the quantity A'^wnlZ-^ov overestimates the real 
WNL/0 ratio by a factor of 1.1 to 7.4, depend- 
ing on the exact age of the burst. As can be seen 
from Meynet (1995, Fig. 6), such a decrease of 
the WNL/0 ratio allows an important steepen- 
ing of the IMF slope, which may bring the values 
closer to a Salpeter-like IMF. 

The possible difficulty of determining the num- 
ber of O stars from rjQ for "evolved" populations 
was already mentioned by V94. One could in 
principle account (to first order) for the evolu- 
tion by simply lowering the upper mass limit to 
the value of the present-day mass function. Such 
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Fig. 1. — Time evolution of r]Q (solid and 
short-dashed line) in an instantaneous burst for 
Z=0.004 and 0.001 respectively. A Salpeter IMF 
from 0.8 to 120 Mq is assumed. To indicate the 
time occurrence of the WR phase the theoreti- 
cal WR/0 number ratio (long-dashed, dotted for 
Z=0.004 and 0.001) is also plotted here. At this 
time the decrease of rjQ with respect to its ZAMS 
value implies an overestimate of the WR/0 ratio 
by typically a factor of ~ 1.1 - 7.4 for the metal- 
licities shown here, depending on the exact age 
of the population 
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a correction has, however, not been applied by 
VC92 and CDC95. To obtain a better quan- 
titative handle on the massive star populations 
in young starburst galaxies the ideal way is to 
synthesize directly the relevant observable quan- 
tities consistently with stellar evolution models. 
Such an approach is presented in the following 
Section. 

2. Evolutionary synthesis models for WR 
galaxies 

We have built population synthesis models 
using the latest Geneva stellar evolution tracks 
(see Meynet et al. 1994 and references therein) 
and atmosphere models from Kurucz (1991) for 
stars with initial masses Mini < 25Mq. To ob- 
tain a reliable description of the ionizing fluxes 
of the hot star population we include the re- 
cent spectra from the "combined stellar structure 
and atmosphere models" ( CoStar) of Schaerer et 
al. (1996ab) and Schaerer (1996) for more mas- 
sive main sequence stars. WR stars are described 
by the theoretical spectra of Schmutz, Leitherer, 
k Gruenwald (1992). Both the CoStar and WR 
atmosphere models account for non-LTE effects 
and stellar winds, which has a strong influence 
on the ionizing flux especially in the He II con- 
tinuum (see Schaerer et al.). 

To study the massive star content in WR 
galaxies we synthesize the following observational 
quantities: 1) H/3 luminosity and H/3 equivalent 
width accounting for stellar and nebular contin- 
uum emission, 2) nebular He II A 4686 emis- 
sion, and 3) broad WR He II A 4686 emission 
feature. The nebular continuum and recombina- 
tion lines were derived assuming N^. = 100 cm~^, 
Te = 10000 K, and solar H/He abundances (see 
e.g. Osterbrock 1989). For the WR feature we 
synthesize the broad He II A 4686 line in order 
to compare our predictions with the high reso- 
lution observations of VC92, which do not in- 
clude broad N III and C III blends at A ~ 46500. 

*At lower resolution these features, often referred to as the 
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Fig. 2. — Predicted He ii/H/3 intensity func- 
tion of the H/3 equivalent width in instantaneous 
bursts at Z =0.008, 0.004, and 0.001. The solid 
lines include both the nebular He ii and the broad 
WR emission. The dashed, dotted and dash- 
dotted lines show the pure WR emission. The 
shaded domain illustrates the region, where a 
strong nebular contribution is predicted. Obser- 
vations from VC92 and CDC95 are identified 
as follows. 1: Mrk 1094, 2: Mrk 712, 3-^4: NGC 
3125 A+B, 5: Tol 35, 6: NGC 1741 B, 7: II Zw 
40, 8: He 2-10 A, 9: Mrk 1236 A, 10: Pox 4, 11: 
Pox 139 
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Consistently with their analysis we adopt an av- 
erage luminosity of 1.7 10^^ ergs~^ for WN stars. 
WC/WO subtypes have a negligible flux in the 
He II A 4686 line. 

Figure |2| shows the predicted evolution of the 
(WR+nebular) He Il/H/3 intensity func- 
tion of the H/? equivalent width for instanta- 
neous bursts assuming a Salpeter IMF ranging 
from 0.8 to 120 Mq. The results for three metal- 
licities Z are plotted. Correspondingly, the ob- 
servations by VC92 and CDC95 are grouped 
into three metalUcity "bins": 0.010 > Z> 0.005 
(open circles), 0.005 > Z > 0.003 (filled trian- 
gles), and 0.003 > Z > 0.0025 (stars). The 
time evolution proceeds towards low 1^(H/?). At 
log VF(H/3) ~ 2.5, the increase of He Il/H/3 corre- 
sponds to the onset of the WR rich phase of the 
staburst. Note the strong decrease of the maxi- 
mum He II/H/3 emission with decreasing Z due 
to strong dependence of the WR population on 
metallicity (cf. Maeder & Meynet 1994). In this 
diagram adopting a flatter IMF would primarily 
translate to an increase of the He Il/H/3 emission 
at a given II^(H/3) by increasing the WR emis- 
sion. 

Figure ^ clearly shows that, except for object 
11, all observations lie within the range covered 
by the model predictions. For each metallicity 
group the agreement is good. Most importantly, 
the observations for objects 1-5 (see figure cap- 
tion) which, according to Meynet (1995) and 
CDC95, would require a very flat IMF [a ~ 1) 
can thus very well be explained with a "stan- 
dard" Salpeter IMF. Such important deviations 
from the non-evolving cluster approach are in- 
deed expected, since according to their Vl^(H/3) 
these bursts have ages i > 4 - 4.5 Myr (see Sec- 
tion |1 . 1[ ) . Our results conflrm the similarly de- 
tailed approach of Kriiger et al. (1992). They 

"WR bump", are hardly separable. Emission in the WR 
bump is larger than in the 4686 feature alone, and the be- 
haviour of its time evolution differs both qualitatively and 
quantitatively from the pure He II 4686 feature modeled 
in this work (cf. Meynet 1995). 



are also in agreement with unpublished observa- 
tions of other WR galaxies (Vacca, private com- 
munication) . 

Only two objects (8 and 11) seem to disagree 
with the Salpeter IMF. He 2-10 (object 8), which 
has a metallicity of Z ~ 0.003, might indeed re- 
quire a somewhat steeper IMF. Pox 139 (object 
11) appears as an extreme object in the VC92 
sample beacause of its large He II/H/? ratio at 
the young age it ~ 2.5 Myr) implied the by high 
value of Ty(H/3). It must, however, be noted that 
VF(H/?) is considerably uncertain due to the very 
low continuum flux measured in the spectrum 
of VC92 (Vacca, private communication). The 
Pox 139 observations of Kunth & Sargent (1983) 
show a He Il/H/3 value lower by a factor of 1.8. 
Better quality observations would be very useful 
for this galaxy. 

2.1. Nebular He ii emission in young star- 
bursts 

The presence of nebular He II emission in gi- 
ant H II regions, requiring a very hard exciting 
spectrum, appears puzzling (see e.g. Garnett et 
al. 1991). The use of the latest CoStar model 
fluxes for O stars and appropriate atmosphere 
models for WR stars (see above) allows us to 
readdress this question here. 

Figure m shows the predicted nebular He Il/H/3 
ratio in instantaneous bursts at different metal- 
licities. First we note that for young bursts we 
predict typical values of /(He Il)//(H/3) between 
5. 10~^ and 2. 10~^. Due to the strong He II con- 
tinuum flux obtained from the O star models, 
which account for non-LTE effects, line blanket- 
ing, and stellar winds (cf. Schaerer et al. 1996b) 
these values are ~ 4 orders of magnitudes larger 
than the values which would be derived from 
other current population synthesis models. 

After ~ 3 Myr the He Il/H/3 ratio increases 
due to the appearance of WR stars. The major 
nebular emission is due to stars in the WC phase. 
Since these stars have a negligible broad stel- 
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Fig. 3. — Predicted nebular emission line ratio 
of He ii/H/3 clS cL function of burst age for differ- 
ent metallicities {Z =0.001 to 0.04) assuming a 
Salpeter IMF from 0.8 to 120 Mq. 

lar emission in the 4686 line, we thus expect an 
important nebular contribution to the total 4686 
emission in the early WC-rich phase. The do- 
main corresponding to this phase is shown by 
the hatched area in Fig. ^. Typical values of 
/(He II)//(H/3) ~ 0.01 - 0.025 are attained. At 
later times (characterized by W^(H/3) < 100 A 
in Fig. the broad stellar component from the 
then important WNL population will dominate 
the 4686 emission. The strongest relative nebular 
contribution is predicted between solar metallic- 
ity and ~ 1/5 Zq , where in the early WR phase 
the nebular flux can make up to ~ 75 % of the 
emission in the He II 4686 line. At Z > 0.02 and 
Z < 0.004 the broad WR feature always domi- 
nates the 4686 emission^. 

At high Z the reduction of the absolute H/3 emission over- 
whelms the increase of the nebular He II/H/3. At low Z, 
due to the low mass loss, the WC population becomes 



First comparisons with observations show the 
following: 1) Two objects in Fig. |2| lie in the do- 
main where the nebular contribution to the to- 
tal WR-l-nebular emission is expected to be high. 
Our prediction is confirmed by inspection of the 
spectrum of Pox 4 from VC92 (their Fig. 6), 
which indeed shows an important nebular con- 
tribution. We expect a similar nebular contri- 
bution in II Zw 40. We note that these objects 
also show the largest electron temperature in the 
VC92 sample, which is an additional confirma- 
tion of their high excitation. For a more firm con- 
clusion about the excitation conditions in Mrk 
712 higher resolution data should be obtained. 
2) Of the 12 objects from Campbell, Terlevich, & 
Melnick (1986) showing nebular He II emission 
four (T 1304-353, T 1304-386, C 1543-^091, F 
30) can well be explained by our models. We sus- 
pect that many of the remaining objects might 
well turn out to show broad He II emission if 
observed with a sufficient signal/noise (see also 
Masegosa, Moles, del Olmo 1991). Their three 
additional objects with broad WR features are 
also compatible with our models. 

3. Summary and outlook 

We have shown that methods, which do not 
account for the evolution of the bursting popula- 
tion will systematically underestimate the num- 
ber of O stars derived from integrated spectra 
for burst ages t > 2 — 3 Myr. As a consequence, 
present WR/0 number ratios in WR galaxies 
(VC92, CDC95) are thus overestimated. There- 
fore IMF determinations based on such number 
ratios will preferentially lead to flat IMF slopes 
(a ~ 1-2; Meynet 1995, CDC95). If one instead 
properly treats the evolution, the observed mas- 
sive star population in WR galaxies is compat- 
ible with resulting from an instantaneous burst 
following a Salpeter IMF, in agreement with pre- 
vious studies using different methods. This con- 
negligible. This explains the dominance of the the broad 
WR emission for these metallicities. 
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elusion may, however, need to be revised if otlier 
systematic uncertainties (photon leakage, dust, 
sht positioning etc.; see e.g. Conti 1993) turned 
out to be important. Obviously it must also be 
noted that the observed properties are probably 
not uniquely described by the case of an instan- 
taneous burst discussed here, and that several 
uncertainties remain (e.g. ionizing fluxes of WR 
stars) . 

Using recent stellar evolution models and new 
spherically expanding non LTE atmosphere mod- 
els for O and WR stars, which in particular yield 
reliable predictions for the ionizing flux in the 
He II continuum, we have proposed a natural ex- 
planation for the observation of nebular He II in 
extragalactic H II regions and alike objects: High 
excitation objects should be intimately linked 
with an important population of WC stars. It is 
indeed feasible to test this prediction by detect- 
ing the broad C IV A 5808 emission line charac- 
teristic of WC stars (see also Meynet 1995). In 
relatively short phases (~ 1 Myr) the expected 
emission (using the WC fluxes from Smith 1991) 
reaches up to /(C IV)//(H/?) ~ 0.02 - 0.2 for 
metallicities between 1/5 Zq and solar. Simi- 
larly, a careful analysis of the C III A 4650 blend 
(cf. Kriiger et al. 1992) and maybe even the 
C III A 5696 line could also be used to measure 
the WC population. 

Although little is presently known about WC 
stars in starbursts the observations of Gonzalez- 
Delgado et al. (1994) showing strong nebular 
He II emission associated with WC stars sup- 
port our picture. Future studies should draw 
their attention on obtaining high quality obser- 
vations (and providing detection/non-detection 
limits !) encompassing both the "traditional" 
WR bump (A 4650 - 4686) with a sufficient spec- 
tral resolution, and the C emission lines. Pre- 
dicted starburst spectra allowing detailed probes 
of WN and WC populations will be presented 
elsewhere (Scha«rer & Vacca 1996, in prepara- 
tion). A systematic study of a large sample of 
young starbursts, including higher metallicity ob- 



jects, will provide crucial tests for stellar evo- 
lution and the understanding of star formation 
processes in these environments. 
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